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Abstract We discuss very low temperature experiments on superconducting micro- 
coolers made of a double Normal metal - Insulator - Superconductor junction. We 
investigate with a high resolution the diflterential conductance of the micro-cooler as 
well as of additional probe junctions. There is an explicit crossover between the single 
quasi-particle current and the phase-coherent Andreev current. We establish a thermal 
model by considering the thermal contribution due to the Andreev current. The related 
increase of the electron temperature is discussed, including the influence of several 
parameters like the phase-coherence length or the tunnel junction transparency. 
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1 Introduction 

The transfer of quasi-particles across the junction between a Normal metal (N) and a 
Superconductor (S) is mainly governed by two processes. 

Single quasi-particles can tunnel from the normal metal to the superconductor if 
their energy E compared to the superconductor Fermi level is larger than the super- 
conducting gap A (E > A). This energy selectivity induces a cooling of the electronic 
population of the normal metal [T] in a S-I-N (where I stands for Insulator) junction bi- 
ased at a voltage below the gap A/e. As the heat current direction does not depend on 
the sign of the bias, S-I-N-I-S micro-coolers based on a double tunnel junction feature 
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a double cooling power and an improved efficiency due to the better thermal isolation 
of the metal. The electronic temperature reduction reaches an optimum at a voltage 
bias just below the gap. In a Al-based device, normal metal electrons can cool from a 
bath temperature of 300 mK down to below 100 mK [T]l2]. 

For an energy E below the gap {E < A), two quasi-particles can tunnel into the 
superconductor and form a Cooper pair in the superconductor. This mechanism is 
called Andreev reflection since it can be viewed as the reflection of an electron into a 
hole [3]- The Andreev current is widely believed to be a dissipation-less current, which 
means that it would contribute only as a charge current. 

The junction normal-state resistance _Rjv is proportional to 1/T, where T is the 
single quasi-particule tunneling probability. In comparison, the Andreev reflection is a 
two-particule process and its probability is proportional to T^. In a ballistic picture [2], 
Andreev reflection is therefore vanishing in a S-I-N tunnel junction, where the interface 
transparency is small. 

In the presence of disorder, the electrons specularly reflected at the superconducting 
interface are confined in the vicinity of the barrier and hit the interface several times. If 
phase coherence is preserved, the probability amplitude of Andreev reflection for every 
attempt add constructively '5' . As a hole and an electron travel on the same trajectory 
but in opposite directions, this addition is immune to the phase randomization induced 
by the disorder. At an energy E compared to the Fermi level in the superconductor, 
the constructive coherent addition of probability builds up over the energy-dependent 
coherence length: 




or by the phase-coherence length L^p if smaller [6]. If one considers a electronic pop- 
ulation at thermal equilibrium and a small bias eV <C ksT, one can define a mean 
coherence length: 

/ hD 

With the effect of phase-coherent confinement taken into account, the Andreev 
channel contributes as a conductance of the order of Rdif f / where Rdiff is the 
resistance of the diffusive phase-coherent normal metal. At low enough temperature 
and bias, it can dominate the conductance of S-I-N junctions with an intermediate 
transparency. A zero-bias conductance peak is then observed, with a width given by 
the Thouless energy of the coherent diffusive region, or the thermal energy 2fc^T if 
larger. The enhanced Andreev current [7l[8] due the phase-coherent confinement by the 
disorder was first observed in Ref. [9]. It was later shown that the Andreev current 
can be modulated by a magnetic flux |10) . In the limit of a strong confinement, for 
example due to a second barrier within the normal metal [11], this effect is more often 
called reflectionless tunneling as the specular reflection at a tunnel barrier appears 
as reduced. 

2 Experimental results 

Fig. [1] left (inset) shows the micrograph of a typical cooler device, where a central 
normal metal Cu island is attached to two superconducting Al reservoirs through tunnel 
junctions. The two 40 nm thick and 1.5 wide superconducting Al electrodes were 
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Fig. 1 (Color online) Left: Differential conductances measured at a cryostat temperature of 
90 mK. Top (red) curve: data of one probe junction 1.55 /im from the cooler junction and 
of normal-state resistance Rpf = 2.76 kl2. The black dotted line is a fit to Eg. [3l describing 
the single quasi-particle tunneling current. Bottom (blue) curve: cooler junction data with a 
normal state resistance Rjv = 1.9 hfi. Inset: Schematics and micrograph of a cooler made of two 
Al-A10(x)-Cu junctions in series. The area of a cooler junction is 1.5 X 0.3 finP . In addition 
to the cooler, one of the three Al-AlO^^-Cu probe junctions on the bottom superconducting 
electrode is visible. The superconducting gap is 2A = 0.43 meV. The voltage axis is normalized 
to A (probe data) or 2 A (cooler data). Right: Normalized differential conductance of the cooler 
junction as a function of voltage at different cryostat temperatures: 240 (purple), 140 (green) 
and 90 mK (blue line). 

in-situ oxidized in 0.2 mbar of oxygen for 3 min before the deposition of the central Cu 
island, which is 4 jj,m long, 0.3 /im wide and 50 nm thick. In addition to these cooler 
junctions, we added three Cu tunnel probes of area 0.3 x 0.3 fim^ on one Al electrode. 
Due to the large volume of the probe Cu electrode, the probe is strongly thermalized 
to the cryostat temperature. In the following, we will describe the experimental results 
obtained on one sample, while we observed a very similar behavior in three more 
samples. At intermediate temperature (above about 200 mK), these samples showed a 
behavior identical to the one reported in Ref. [T3] . The charge current in both junctions 
can then be described by the sole single quasi-particles contribution. The cooler shows 
an electronic cooling, while the probe does not, as expected. 

Let us now concentrate on the very low temperature regime. Fig. [l]left shows the 
differential conductance of the cooler double junction and of one of the probe junction 
at the cryostat temperature of 90 mK. It was obtained by numerical differenciation of 
the measured current-voltage characteristics. At intermediate bias, the probe junction 
data is well fitted (black dotted line) by considering only the single quasi-particle tunnel 
current given by: 



1 f 

lTiV) = -^ / Ns{E){fN{E~eV)^fN{E + eV)]dE, 




(3) 



[14j with an electronic temperature of 105 mK, which is slightly higher than the cryostat 
temperature. 
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Close to zero bias, the cooler and the probe feature a peak in the differential 
conductance. This effect cannot be accounted for by a linear leakage as it would lead 
to a saturation of the conductance near zero bias. As discussed in Ref. [14], this zero 
bias anomaly cannot be fitted by considering a non-equilibrium distribution in the 
normal metal or by considering a smeared density of states in the superconducting 
electrodes |15] . Below about 200 mK, the zero bias conductance increases when the 
cryostat temperature is lowered (see Fig. [U right), which suggests that it is a phase- 
coherent effect. In the cooler junction, the differential conductance is decreased at 
intermediate bias and more peaked at the gap edge compared to the one of the probe 
junction, which exemplifies the cooling of the normal metal island electrons. 



3 The Andreev current at thermal equilibrium 

In the following, we will ascribe the zero bias enhancement of the differential conduc- 
tance to a phase-coherent Andreev current. We will use the theory of Ref. [8lll6]. In 
our samples, the coherence length Lx at T = 90 mK is about 0.33 fim, which is of the 
order of the junction dimensions. At low energy, the propagation of the Cooperon in 
the electrode is cut by the phase-coherence length Ltp, which is expected to be about 2 
fim, i.e. larger than the junction dimensions. Therefore, we will use the ID regime for 
the diffusion of the Cooperon in both the normal metal and the superconductor elec- 
trode. We take into account the finite gap of superconductor, so that the calculation 
is valid for eV, kT < A, and we include the disorder both in the normal metal and in 
the superconductor [17] . 

We have fitted the probe data by calculating the sum of the single quasi-particle 
current It and the phase-coherent Andreev current 7^. Fig. [2] left shows the excellent 
fit with the two contributions of the Andreev current at low bias and of the single 
quasi-particle current at higher bias. The fit parameters are: L^p = 1.5 fim, A = 0.228 
meV, T = 105 mK. We took the measured values of the diffusion coefficient D = 80 
cm^/s and R]\[ = 2.76 kJ7. We have used the same parameters to fit successfully the 
experimental data from the two other probe junctions. In the fit, we had to scale the 
phase-coherent Andreev current by a multiplying factor AI = 1.37. This factor [TO] 
could be due to small inhomogeneities in the tunnel barrier, which are not considered 
here. 



4 The behavior of the cooler junction in the presence of an Andreev 
current 

In order to understand the behavior of the cooler junction, we need to consider the heat 
balance in the central normal metal. Here we assume a quasi-equilibrium situation: the 
electrons and the phonons in the metallic island follow a thermal distribution function 
at a respective temperature Te and Tp/j, which are in general different from the bath 
temperature Tftaj^ of the cryostat. 

The single quasi-particle current is responsible for the cooling power: 

Pcool = ^)ns{E)[fN{E - ^) - fsiE)]dE (4) 
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Fig. 2 (Color online) Experimental current- voltage characteristics (full red lines) of the probe 
(left) and cooler (right) junctions at a cryostat temperature of 90 mK, compared to calculated 
curves. The dotted lines are fits with the calculated phase-coherent Andreev current and the 
dashed lines shows the fit with the calculated single quasi-particle current. The parameters 
are: D = 80 cm^/s, = 1.5 fim, M = 1.37 (probe) or 0.49 (cooler), 2A = 0.43 meV, KA = 
144 W.K-4. 



The cooling power is compensated by the electron-phonon coupling power 

Pe-ph = SU{t! -T^h), (5) 

so that 2Pcool + Pe-ph ~ 0, where the factor 2 is due to the fact that we have N-I-S 
junction in series. We consider that the heat given to the normal metal phonons is 
compensated by the Kapitza coupling with the phonons of the substrate kept at the 
bath temperature T^^ath- Pr ~ ^^iTj^ath ~ '^ph)^ ih&t P^^ph + Pk = 0. Here P^ is 
the power flow through the Kapitza resistance, K is the Kapitza parameter depending 
on the materials in contact, and A is the contact area. The Kapitza thermal resistance 
is significant for intermediate temperatures and above (T > 300 mK), which can lead 
to the cooling of the normal metal phonons [13]. We have assumed that the Andreev 
current does not dissipate any heat. 

To fit the experiment, we first solve numerically the thermal model discussed above 
for a relatively high temperature (T > 300 mK), where the contribution of the phase- 
coherent Andreev current is negligible. We take the electron-phonon coupling coefficient 
S = 2 nW./im~^.K~^ and obtain from the fit the Kapitza coupling parameter KA = 
144 W.K~^. The Kapitza coefficient K is comparable to the one found in our previous 
experiments [IS]. 

We turn afterwards to the very low temperature regime of interest here. Fig. fright 
shows the direct current-voltage characteristic obtained from the cooler junction (full 
red line) along with calculated curves (dashed and dotted lines) . The dashed line is the 
calculated current-voltage characteristic at a 90 mK cryostat temperature including 
the charge and heat currents of the single quasi-particle tunneling only. The agreement 
is poor, which confirms the need to include the Andreev current contribution to the 
junction. The dotted line shows the result of the calculation based on the thermal 
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Fig. 3 (Color online) Current-voltage characteristic of the cooler junction (full red line) along 
with curves calculated with the calculated thermal model and including an excess leakage due 
to a linear resistance of 6 (dash-dotted line), 10 (dashed line) and 100 MJ7 (dotted line). 

model with the charge current given by the sum of the single quasi-particle current 
and the Andreev current. The fit parameters are the same than for the probe except 
for the scaling factor M = 0.49. The difference with the probe junction factor is not 
understood, although it could be due to the difference in geometry between the two 
junctions. The addition of the phase-coherent current provides an acceptable fit at low 
bias but shows a clear discrepancy at intermediate voltage. The experimental curve 
predicts a larger current than what is obtained from the thermal model. This demon- 
strates that an excess dissipation term or an extra current contribution is missing in 
the thermal model. 

As a possible explanation, Fig.|3]shows the comparison between the experiment and 
a calculated curve from the thermal model, which includes an additional dissipation 
due to a linear resistance. Here the leakage contributes both as a dissipation in the 
normal metal and as a current across the junction. The differential conductance at 
zero bias of the cooler junction (see Fig. [l| gives a minimum leakage resistance of 20 
MJ7. Fig. 13] shows that adding such an extra dissipation term to the thermal model 
discussed above does not provide a good description of the experiment. 



5 The Andreev current induced dissipation 

Let us now discuss the heat transfer due to the Andreev current. The work performed 
by the current source feeding the circuit with the extra Andreev current Ij\ generates a 
Joule heat ~ IaV that is deposited in the normal metal [18] . This heat is deposited 
entirely in the normal metal and does not perturb the superconductor. Hence, the net 
cooling power of the S-TN junction is reduced and can be re-defined as Pcool ~ Pa- 

Fig.|4]shows the quantitative comparison at 100 mK of the cooling power Pcool due 
to single quasi-particle tunneling (red line) and the dissipation P4 due to the Andreev 
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Fig. 4 (Color online) Calculated power due to quasi-particle cooling Pcool (red line) and the 
heat dissipation due to the Andreev current (blue line) as a function of voltage bias for a 
N-I-S junction at T^^ath = 100 mK. The parameters used in the calculation are the same as in 
Fig. [2] 

current (blue line). Near zero bias, Pcool is almost zero due to the absence of quasi- 
particles. It attains its maximum near the gap. The Andreev current induced dissipation 
Pa increases sharply near zero bias. Close to the gap voltage, the cooling power out- 
does the Andreev dissipation. As the latter depends strongly on the transparency of 
the junction, it surpasses the single quasi-particle cooling at a varying temperature 
and bias. In the present device with a tunnel barrier transparency of about 10~^, the 
Andreev current dissipation becomes relevant only below 200 mK. 

Fig. [5] shows the schematic of the full thermal model of the device. Here we have 
included the work Ij^V done on the central metallic island by the current source. At 
steady state, the heat balance for the normal metal electrons can be rewritten as: 

"iPcool + Pe-ph - IaV = 0. (6) 

With this complete heat balance equation taken into account, we solve the thermal 
model and calculate the total contribution to the current in the cooler junction. Fig. |6] 
left shows the comparison of the experiment (full colored lines) and the thermal model 
(dotted lines). The agreement is very good at every accessible cryostat temperature. 
The Andreev current contribution becomes negligible compared to the single quasi- 
particle current contribution at temperatures above 200 mK. 

The above conclusion on the Andreev heat is independent of the phonon cooling 
[13) . Assuming the perfect thermalization of the phonons to the substrate temperature 
would change the total calculated current by less than 2 % at 90 mK, which means 
that phonon cooling has a negligible role in the data analysis at very low temperature. 
This is consistent with the expected negligible amplitude of the phonon cooling in this 
temperature range. 

The above fit with the thermal model also provides us with the electron temper- 
ature for every bias. Fig. [6] right shows the calculated electron temperature in the 
central metallic island as a function of voltage bias across the cooler junction and for 
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Fig. 5 Schematic of the complete thermal model of S-I-N-I-S micro-coolers at very low tem- 
perature (see text). 
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Fig. 6 (Color online) Left: Current voltage characteristic of the cooler junction at different 
cryostat temperatures T^^jtij together with the calculated best fit from the full thermal model 
including the charge and heat contribution due to the Andreev current. The temperatures 
are from top to bottom 430 (purple), 330 (orange), 230 (green) and 90 mK (red line). Right: 
Dependence of the calculated electronic temperature with the voltage with the parameters 
obtained from the fit to the experiment and for a series of cryostat temperatures: 230 (green 
line), 140 (blue) and 90 mK (red). 



different cryostat temperatures Tfjath = 90, 140 and 230 mK. At very low temperature, 
the electron temperature first increases with the bias as the Andreev current- induced 
dissipation is dominant. When the voltage bias approaches the gap, the single quasi- 
particle tunneling based cooling dominates. As the bath temperature increases, the 
Andreev current induced dissipation becomes less efi'ective and for Tf,ath ~ 230 mK 
the electronic cooling always prevails [19| . 

Although the Andreev reflection is a small effect as a charge current in S-I-N- 
I-S cooling devices, the related heat contribution is extremely efficient at very low 
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temperature. A basic explanation for this is the following. At very low temperature the 
quasi-particle based cooling power has a small efficiency compared to the Joule power 
IV. It is of the order of Te/A, which is about 5% at a 100 mK electron temperature. 
In contrast, the Andreev current induced dissipation is the full Joule power. 

6 Conclusion 

As a summary, we observed a peak in the differential conductance at low bias in the 
probe and cooler junctions of superconducting micro-coolers, which is due to a phase- 
coherent Andreev current. A quantitative thermal analysis of the micro-cooler behavior 
demonstrates the importance of the dissipation induced by the Andreev current, as it 
dominates the cooling power at very low temperature over a significant voltage range. 

The above conclusion poses a challenge to diminish the dissipation induced by the 
phase-coherent Andreev current in S-I-N-I-S micro-cooling devices. What are the right 
parameters needed to optimize the cooling power? The induced dissipation due to 
the Andreev current depends strongly on the transparency. It scales as 1/ Rj^ whereas 
the quasi-particle cooling Pcool scales as 1/ Rn ■ A resistance optimum in terms of cool- 
ing power can then be found for every temperature. As the temperature decreases, 
the high-energy tail of the electron energy distribution is vanishing, which reduces the 
cooling. In contrast, the Andreev current would increase. It is thus expected that the 
optimum resistance increases when the electronic temperature decreases. This qualita- 
tive discussion needs to be completed by a more quantitative calculation. 

For a fixed tunnel resistance, a smaller phcise-coherence length Lip would diminish 
the phase-coherent Andreev current as well as the related dissipation. For instance, a 
disordered material would bo a better choice as a normal metal since it would have 
a much shorter phase-coherence length than a pure metal like Copper. The latter 
strategy seems clearly promising for new devices with an improved efficiency at very 
low temperatures. 
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